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Abstract 

The near-field interaction between fluorescent emitters and graphene exhibits rich 
physics associated with local dipole-induced electromagnetic fields that are strongly 
enhanced due to the unique properties of graphene. Here, we measure emitter lifetimes 
as a function of emitter-graphene distance d, and find agreement with a universal 
scaling law, governed by the fine-structure constant. The observed energy transfer- 
rate is in agreement with a 1 /d 4 dependence that is characteristic of 2D lossy media. 
The emitter decay rate is enhanced 90 times (transfer efficiency of ~ 99%) with respect 
to the decay in vacuum at distances d ~ 5 nm. This high energy-transfer rate is mainly 
due to the two-dimensionality and gapless character of the monoatomic carbon layer. 
Graphene is thus shown to be an extraordinary energy sink, holding great potential for 
photodetection, energy harvesting, and nanophotonics. 

*To whom correspondence should be addressed 

tlCFO 

*CSIC 

1 



Graphene, a genuinely two dimensional material composed of a single layer of carbon 
atoms, has rapidly generated great interest since its experimental isolation^ thanks to its ex- 
traordinary electronic, optical, and mechanical properties. Due to its gapless band structure 
and linear electronic energy dispersion,^ graphene exhibits frequency-independant light ab- 
sorption over a broad spectral region in the visible and infrared,^ which is governed only by 
fundamental material-independent constants: its absorbance is given by na ~ 2.3%, where 
a ~ 1/137 is the fine-structure constant.^ This new material displays high room-temperature 
mobilities^ up to 250,000 cm 2 V _1 s _1 , electrically tunable carrier concentration, and bipolar 
field response.^ These optical and electronic properties have promoted the use of graphene 
for a multitude of opto-electronic applications 6 and as a potential platform to exploit light- 
matter interactions under ambient conditions.^ 

Recently, efficient energy-transfer from bio-molecules to graphene has sparked tremen- 
dous interest for sensing purposes.^^ From a fundamental point of view, the near-field 
interaction between an emitter and a purely two-dimensional material is particularly inter- 
esting because it allows for the exploration of new limits of light-matter interactions.^^ 
First of all, due to the two-dimensional and gapless character of graphene the magnitude 
of the non-radiative coupling is strongly enhanced relative to other lossy materials. This 
strong coupling has been predicted to produce substantial energy transfer up to distances of 
30 nmP Additionally, the non-radiative energy transfer rate T nr as a function of distance d 
between the emitter and graphene has been predicted 16 * 17 ! to scale as <i -4 , in contrast to the 
energy transfer of emitters near conventional semi-infinite materials, where the magnitude 
scales with g? -3 . Finally, the emitter decay rate enhancement follows a universal scaling law 
governed by the fine-structure constant and the ratio of d and the emitter wavelength A.^ 

From an application point of view, we envision novel types of hybrid systems that combine 
graphene with strong light absorbers/emitters (e.g., quantum dots and fluorescent molecules) 
to enhance the intrinsically weak graphene absorption, and thereby, also the device efficiency 
for light harvesting and photodetection. Finally, the scaling of the decay rate as a function of 
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distance makes the emitter-graphene system promising as a nanoscale ruler, as proposed in 
Ref. 17. Further advances towards these new effects and applications require a quantitative 
experimental study, which has not yet been realized. 




Rhodamine Graphene 



Figure 1: (A) Energy diagrams of an optical emitter and a nearby graphene monolayer 
sheet. Upon optical excitation of the emitter, the relaxation of the excited state occurs 
for short distances primarily through non-radiative decay by dipolar coupling to electron- 
hole pair transitions in the carbon layer, and to a lower extent through the emission of 
radiation. This leads to emission quenching and a shorter lifetime of the emitters. (B) 
Schematic representation of the sample structure. The separation between the emitters and 
the graphene flake is implemented through iterative atomic layer deposition of a TiC>2 spacer 
layer, as explained in the text. 

The efficient energy transfer between light absorbers/emitters and graphene relies on 
the strength of their near-field interaction. Recent theoretical studies^^ suggest that this 
interaction is mediated by non-radiative coupling between the emitter dipole and electron- 
hole pair excitations in graphene (i.e., Forster-like energy transfer, see Figure which in 
turn results in higher decay rates and emission quenching (i.e., the energy released from the 
emitter, that would otherwise be re-emitted as light, is absorbed by the graphene). This 
picture is consistent with recent experiments that have demonstrated that the fluorescence 
of nitrogen vacancy centers in diamond,^ CdSe/ZnS nanocrystal emitters,^ fluorescent 
,~[y es EEE2] anc [ dye-labeled DNA^l is quenched by the presence of graphene. 

In this Letter, we perform a quantitative study of the near-field interaction between 
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graphene and nearby emitters. Specifically, we measure the decay rate of a thin layer of 
rhodamine fluorescent molecules (emitters) coupled to a monolayer graphene flake as a func- 
tion of emitter-graphene distance. We show that graphene induces a dramatic change in 
emitter lifetime, increasing the decay rate of the excited molecules by up to 90 times. This 
corresponds to >99% of the energy stored in the emitters being transferred to graphene. We 
compare the results with a simple (but rigorous) analytical model and find good agreement 
without any fitting parameters. The observations reveal that the strong near-field inter- 
action originates from the unique properties of graphene: its gapless and two-dimensional 
character as well as its charge carriers being relativistic massless Dirac fermions. As a result, 
the non-radiative energy transfer to graphene is governed by the fine-structure constant and 
exhibits a 1/ <i~ 4 universal scaling with the distance d between emitters and graphene. 



Experiment 

We measure the optical-excitation lifetime of a thin layer of rhodamine molecules (< lnm 
thick 24 ) by probing spatially and temporally resolved fluorescence using a home-built confo- 
cal microscope. A pulsed green laser (532 nm wavelength) at 40 MHz excites the molecules 
and their fluorescence is recorded with an avalanche photodiode detector (APD). Through 
time-correlated photon counting with a Picoharp time correlator, we obtain the lifetime of the 
emitters. The sample consists of a monolayer graphene flake placed on top of a Si-SiC>2 (285 
nm) substrate. The graphene is covered by a TiC>2 spacer layer, realized through atomic layer 
deposition, on top of which we deposit the rhodamine layer, followed by a final capping layer 
of PMMA that provides stability of the emitters (see Figure [TJ3). We modify the distance 
between the emitters and the graphene by increasingly depositing additional layers of TiCV 
In this way, we measure the lifetime of the emitters in the vicinity of a single graphene flake, 
while gradually increasing the graphene-emitter distance (i.e., the TiC>2 layer thickness). We 
choose experimental parameters such that processes other than non-radiative coupling to 
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graphene are much less likely to occur. First, the spacer layer thickness (5 — 20 nm) is much 
smaller than the emission wavelength of the rhodamine molecules (650 nm) and is electri- 
cally isolating, so that charge transfer between the emitters and graphene is prevented.^ 
Second, we use graphene with low intrinsic doping such that plasmons are not excited at this 
wavelength.^ A more detailed description of the experimental setup and sample fabrication 
is presented in the Supporting Information. 
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Figure 2: (A) Optical microscope image of the single layer graphene flake (confirmed by Ra- 
man measurements) used in the experiments. (B) Fluorescence image of rhodamine molecules 
in the region depicted in panel (A). The distance separating between the emitters and the 
graphene flake is 11 nm. Fluorescence quenching from the emitters on top of graphene is 
clearly observed. (C) Two emitter lifetime curves on a region above graphene (red) and not 
above graphene (blue). The solid lines are bi-exponential fits to the data. (D) Lifetime image 
of the region depicted in panel (A), showing a clear reduction of the rhodamine lifetime in 
the region above the graphene flake. 
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Results and discussion 



Figure [2|A. shows an optical microscope image of the graphene flake used in the experiments, 
which can be correlated with Figure [2j3, containing a fluorescence image of the same area. 
Clear fluorescence quenching of the emitters that are located 11 nm above the graphene sheet 
is observed, in agreement with earlier experimental works.^H22 This fluorescence quenching 
is due to the non-radiative energy transfer processes from the emitter to graphene. 

Measurements of the lifetimes of the emitters respectively on top of graphene and outside 
graphene (above the substrate) are shown in Figure [2p. As expected, the lifetime for emitters 
on graphene (defined as T) is shorter compared with emitters on the substrate (defined as 
r s ). Quantitative analysis (see Supporting Information) shows that there are two dominant 
contributors to the emission: the rhodamine emitters and a small background from the 
PMMA capping layer. Therefore, we fit the lifetime data with a double exponential decay 
A e -t/m ho _|_2?e~ t//rpMMA , where t is the delay time between laser excitation and fluorescence, 
TRho and tpmma are the lifetimes of the rhodamine emitters and the PMMA protective layer, 
respectively, and A and B are the corresponding fluorescence intensities (in count rates). We 
find that, for all graphene-emitter distances, the emission from rhodamine is higher than 
that of the PMMA layer (A > B) and that, due to the large thickness of the PMMA layer 
and therefore smaller coupling between graphene and PMMA, tpmma is independent of the 
distance, yielding a constant lifetime of tpmma = 1.8 ns. The solid lines in Figure |2p are fits 
according to this double exponential decay and yield rhodamine decay rates r s = (3.5ns) -1 
on the substrate and Y = (1.0ns) -1 on graphene, which give a decay rate ratio T /T s of 3.5 for 
a graphene-emitter distance of 11 nm. In Figure [2p we show the extracted emitter lifetime 
r Rho for the same region as the optical and fluorescence images of Figure [2]A.,B. The good 
correlation between all three images confirms that it is the graphene flake that causes the 
enhancement of the decay rate of the emitters. 

In order to produce a complete quantitative analysis and show the remarkably strong 
coupling between graphene and the emitters, we investigate how the ratio r/T s depends 
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Figure 3: Emitter lifetime reduction as a function of distance between the emitters and 
graphene. (A) Lifetime curves for emitters on top of graphene for distances of 5nm, 10 nm, 
and 11 nm, as well as for emitters outside the graphene. Solid lines are bi-exponential fits to 
the data. (B) Dependence of the rhodamine lifetime as a function of the distance separating 
them from the graphene. (C) Decay rate enhancement T/r s obtained from the lifetime 
measurements as a function of distance, yielding up to a factor 90. Dashed lines: analytical 
simulation for point-dipole emitters with their orientation parallel (Tii) and perpendicular 
(r_|_) to the graphene plane. The red solid line represents the weighted average over dipole 
orientations, (l/3)r_|_ + (2 /3) T || . The green line represents the calculated r/T s (weighted 
average over dipole orientations) for the multilayer structure which includes the Si (0.5 mm), 
SiC>2(285 nm), TiC>2 (variable thickness), and PMMA (50 nm), with dielectric constants 
of 14.8, 2.12, 5.7, and 2.22 respectively.^ (D) Energy transfer efficiency as a function of 
distance, calculated based on the lifetime measurements presented in (C). 
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on the emitter-graphene distance (see Figure [3]). In panel A, we present three lifetime 
measurements with a spacer layer of 5nm, 10 nm, and 11 nm, respectively, as well as the 
lifetime measured on a region without graphene, showing a clear reduction in lifetime as 
the emitters are placed closer to the flake. We repeated the measurement of the emitter 
lifetime outside the graphene flake region for all thicknesses of the spacer layer d and found 
a lifetime of 3.5 ± 0.3 ns that did not depend on the thickness of the Ti02 spacer layer. The 
analysis indicates that the lifetime of the emitters increases with increasing distance until 
the emitters and graphene are separated by 20 nm (Figure [3^3), where the measured lifetime 
contrast between emitters on graphene and outside graphene disappears. The corresponding 
ratio r/r s , presented in Figure |3p, shows remarkably strong lifetime modifications reaching 
up to 90 for a distance of 5 nm, which is an extraordinary result considering that this lifetime 
reduction is due to a single layer of atoms. We calculate the energy transfer efficiency via 
the relation rj = 1 — r s /T, as is customary in Forster-like processes. We find very efficient 
energy transfer, yielding more than 99% at short distances, as depicted in Figure |3]D. We 
also note that the high efficiency of the energy transfer spans over a large distance range, 
with values > 85% at distances up to more than 10 nm between the emitters and graphene. 



Comparison with theory 

We compare our experimental results for the distance scaling of the energy transfer with a 
semi-classical model of an emitter (energy donor) coupled to a neighboring material (energy 
acceptor). The emitters are described as classical dipoles and the emission rate is worked 
out from the power that these dipoles radiate.^ Details of the model are given in the 
Methods. For tutorial purposes we first consider an emitter placed in vacuum at a distance 
d from graphene supported on a medium of permittivity e, while rigorous results for the 
more complicated multilayer structured of the experiment are given below. In the range 
of distances d — 1 — 15 nm, we find the decay rate T g from emitter to graphene to be well 
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described by (see Methods) 

where Ao is the free-space emission wavelength, To is the rate in vacuum, and v — 1 {y = 2) 
when the emitter dipole is oriented parallel (perpendicular) to the graphene. Interestingly, 
the energy transfer rate does not depend on any material-dependent parameter related to 
the graphene.^ The origin of this universal scaling of the energy transfer between an optical 
emitter and graphene lies in the universal value for the optical conductivity e 2 /h, which is 
attributed to the gapless and 2D character of the lossy graphene system, as discussed in the 
Methods section. 

In Figure |3p we compare this theoretical result with the experimentally obtained lifetime 
data. In order to yield a more realistic comparison, we have straightforwardly extended the 
theory to include the multilayer environment of our samples as well as retardation effects 
following the methods of Ref. 27. We find that the calculated lifetime ratio r/r s from this 
extended theory (green line in Figure 3c) is comparable to Eq. ?? for d > 4 nm (red line in 
Figure 3c). We find that theory and experiment are in excellent agreement both qualitatively 
and quantitatively without any fitting parameters. 

Both the measured data and Eq. ?? show the expected d~^ dependence for the energy 
transfer rate to a two-dimensional material. It is well established that Forster-like energy 
transfer has a typical scaling law d~~ n , where n is determined by the dimensionality of the 
system.^ For two single-point coupled dipoles, the emission rate follows a d~ Q dependence. 
If one of the dipoles is replaced by a line of dipoles, integration reveals a d~ 5 scaling. Anal- 
ogously, for a two-dimensional array of dipoles, a <i -4 scaling is obtained, and finally Forster 
processes scale with <i~ 3 when a single dipole interacts with a bulk of dipoles. Therefore our 
result is in agreement within the error with the expected scaling, given the two-dimensionality 
of the graphene sheet. 
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The observed energy transfer rate of an emitter close to graphene is very strong and 
reaches up to 90 times the decay rate in vacuum for a distance of 5nm. This can be 
attributed to two factors: the two-dimensional character of graphene and the relatively high 
value for its optical conductivity (strictly speaking the real part of its conductivity, which 
accounts for ohmic losses). The latter is in particular high for graphene due to the absence 
of a band gap which leads to a stronger coupling between electron-hole pairs and optical 
fields, compared with materials which have a bandgap. This high decay rate in front of an 
atomically thin material is counterintuitive, as one would expect that bulk, lossy materials 
should be more effective in absorbing from a nearby emitter. 

To put these results for graphene in context, we compare the decay rates of an emitter 
coupled to graphene and other materials which have a semi-infinite or (hypothetical) thin- 
film configuration. The model for the other materials starts from the same general form as 
for graphene (see Methods section). For a semi-infinite material with permittivity e m , we 
have (with the emitter in vacuum) 

r 3D /r ~ i + ^iml^jf^V, (2) 

6D/ 128vr 3 \e m + l}\ d J ' V ; 

while for thin films, we find 



r 2D /r ~ l + ^^Im(e m -l/e m )t^, (3) 

where t and e m are respectively the thickness and permittivity of the film . These equations 
reveal the <i -3 and <i -4 scaling for bulk and 2D-materials. 

Tutorially, and n order to provide a more quantitative comparison, we calculate the decay 
rates of emitters close to both graphene, semi-infinite materials and thin films of 1 nm. Figure 
[4] shows these decay rates for three different materials: graphene, silicon, and InSb. The 
latter is a highly absorbing semiconductor because of its small band gap. These calculations 
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Figure 4: Decay rate enhancement versus distance for graphene, low-bandgap bulk InSb, 
and thin films of InSb and Si. The calculations are based on Equations (1) for graphene and 
Equation (2) for bulk materials, while for the thin films, we use a similar procedure as in 
the multilayer system.^ 
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provide us with two important insights. First, for smaller distances (< 10 — 50 nm, depending 
on the material), thin films lead to higher energy transfer rates than bulk materials. This 
is due to the crossover of the <i -4 scaling for 2D systems compared with the <i -3 scaling 
of bulk materials. Second, a (hypothetical) thin-film of a low band gap material shows 
similar energy transfer rates as graphene. This is because both graphene and low band 
gap semiconductors are very strong light absorbers and exhibit similar values for the real 
part of their optical conductivities. From the perspective of applications, our results show 
the interesting prospective of enhancing the energy transfer to a material by reducing its 
dimensionality, which can be readily achieved by using graphene. We note that for increasing 
distances, the optimal material thickness for highest energy transfer is a slightly thicker film 
than a monolayer, with a linear dependence of this optimal thickness on distance. 

Conclusion 

Our results show a strong enhancement of the energy transfer rate of emitters coupled to 
graphene as they come closer to the carbon layer, with a decay rate enhancement factor 
of up to 90, and energy transfer efficiencies above 85% for a distance up to 10 nm. The 
experimental results are quantitatively consistent with a rigorous model based upon non- 
radiative energy transfer revealing a universal scaling of the decay rate of an emitter coupled 
to graphene, governed by the fine-structure constant a and the ratio d/X. 

The coupling between fluorescent materials and graphene is remarkably strong due to its 
two-dimensionality and gapless character. This can be exploited in numerous applications of 
varied nature. Noting that the dependence of the decay rate of the emitters versus distance 
is governed by material-independent fundamental constants and substrate optical properties, 
the graphene-emitter fluorescence can be used as an absolute nanoscale ruler Taking into 
account graphene's high mobility, it is foreseeable that combining highly absorbing emitters 
with graphene will result in highly efficient photo-detection and energy harvesting devices. 
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In photonics, light can be shed on the dynamics of dark molecules (i.e., molecules with low 
quantum efficiency, whose ratio of intrinsic non-radiative to radiative decay rates is large). 
For example, if a dark molecules is coupled to graphene with an energy transfer rate that is 
larger than the intrinsic non-radiative decay rate, energy from the excited state of the dark 
molecule can be extracted, allowing for the study of excited state dynamics of dark molecules. 
Finally, the measured long range of the energy transfer can be used through DNA-length 
dependent fluorescence for biomolecule sensing applications with nanometer resolution.^ 



Methods 

We use a semi-classical model of emitters coupled to a nearby material. The emitters are 
described as classical dipoles and the emission rate is worked out from the power that these 
dipoles radiate.^ This produces results that are in agreement with a quantum-optics anal- 
ysis.^ For simplicity, we consider in this section the emission from an emitter placed in 
vacuum above an absorbing material. However, we have carried out a straightforward exten- 
sion of this analysis to calculate the decay rate of an emitter in a multilayer system as that 
of the experiment (see Figure The decay rate T for an excited emitter in front of an 
acceptor material in the long-wavelength limit can be expressed as an integral over parallel 
wave vectors ku asP 

r/r ° = 1 + m r e ~ 2fc " dim ^}> ^ 

where To is the rate in free space, r p is the k» -dependent Fresnel reflection coefficient for 
p-polarized light, Ao is the light wavelength, and v — 1 (v — 2) when the dipole is oriented 
parallel (perpendicular) to the surface. The distance dependence of the decay rate is related 
to the fen dependence of r p . 

For graphene, we have Im{r p } = Im{ e+1+ 4^fc /^ }) where a is the graphene conductivity, 
to is the photon frequency, and e is the substrate permittivity. One can approximate^ a = 
e 2 /4:h, so that Eq. ?? admits a closed-form analytical solution T/Tq = 1 + uC I(x), where C = 
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3(e+ l) 2 /2(na) 3 is a constant and x = [4(e + l)/a](<2/ Ao). Here, I(x) = l/x 2 + Ci(x) cos(x) + 
si(x)sin(x) is a function related to the Ci and si functions (tabulated in Ref.^) that can 
be approximated (within < 5% maximum relative error) as I(x) ~ l/(x 2 + x 3 /3 + x 4 /6). In 
the 3 — 15 nm range, this yields the result of Eq. ??, where we find a T ~ dependence 
(similar to a conventional thin film). At small distances d < 3aAo/[4(e+ 1)] ~ lnm (below 
the range of our measurements), the first term in the denominator of I(x) dominates, giving 
rise to a T ~ dr 2 dependence. At even smaller distances d < vp/u ~ 0.3 nm, where the fcii 
dependence of a becomes relevant (nonlocal effects), one recovers a ~ <i~ 3 behavior typical 
of a semi-infinite medium. 

For other materials, in the emission wavelength under consideration and at distances be- 
low ~ 20nm, small k\\ components dominate the integral of Eq. ??, and we can approximate 
r p m gkl} (n = 0, 1, 2, or 3, depending on the dimensionality of the system), which leads to 

3i/Ag (n + 2)!Im{g} 
r/r ° = 1+ 32^ (2d)«+3 • ( 5 ) 

Thus, the dominant distance-dependent term depends on the type of system under consid- 
eration. We examine the 3D and the 2D case. We find that in the 3D case (n = 0) the 
Fresnel coefficient behaves as r p — g — (e — l)/(e + l), where e s is the permittivity of the 
substrate, leading to T /To = 1 + C^D(^o/d) 3 , where C^d = 3^g/1287r 3 . In contrast, for a 
thin film of permittivity e and thickness t <C d, we find r p = [(e 2 — l)/2e]fciit (i.e., n = 1), 
leading to T/Tq = 1 + C2DtA()/<i 4 , where C2D — (9z//5127r 3 )Im(e— 1/e). The latter result has 
the same distance dependence as we find in graphene, which can then be ascribed to the 2D 
character of the carbon layer. In fact, we find from the above expression for r p in graphene 
that we can approximate Im{rp} w (2irak\\/ko)/(e + l) 2 , where ko = ui/c is the free-space 
light wave vector, under the assumption that irak^/ko <e+l, thus recovering the same 
n = dependence as for an absorbing thin film. 

The theory provided in Figure [3] is obtained without any fitting parameters. We find 
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good agreement by representing the graphene through its DC conductivity a = e 2 /Ah. 

Acknowledgment. This work has been supported in part by the Fundacicio Cellex 
Barcelona, the ERC Starting grant no. 307806 (CarbonLight), the ERC Career integration 
grant 294056 (GRANOP), the Spanish MICINN (MAT2010-14885 and Consolider Nano- 
Light.es), and the European Commission (FP7-ICT-2009-4-248909-LIMA and FP7-ICT- 
2009-4-248855- N4E). L.G. acknowledges financial support from Marie-Curie International 
Fellowship COFUND and ICFOnest program. K.J.T. acknowledges financial support from 
NWO Rubicon grant. 

References 

(1) Novoselov, K.; Geim, A.; Morozov, S.; Jiang, D.; Zhang, Y.; Dubonos, S.; Grigorieva, I.; 
Firsov, A. Science 2004, 306, 666-669. 

(2) Neto, A. H. C; Guinea, F.; Peres, N. M. R.; Novoselov, K. S.; Geim, A. K. Rev. Mod. 
Phys. 2009, 81, 109-162. 

(3) Mak, K. F.; Sfeir, M. Y.; Wu, Y.; Lui, C. H.; Misewich, J. A.; Heinz, T. F. Phys. Rev. 
Lett. 2008, 101, 196405. 

(4) Nair, R.; Blake, P.; Grigorenko, A.; Novoselov, K.; Booth, T.; Stauber, T.; Peres, N.; 
Geim, A. Science 2008, 320, 1308. 

(5) Mayorov, A. S.; Gorbachev, R. V.; Morozov, S. V.; Britnell, L.; Jalil, R.; Pono- 
marenko, L. A.; Blake, P.; Novoselov, K. S.; Watanabe, K.; Taniguchi, T.; Geim, A. K. 
Nano Lett 2011, 11, 2396-2399. 

(6) Bonaccorso, F.; Sun, Z.; Hasan, T.; Ferrari, A. Nat. Photon. 2010, 4, 611-622. 

(7) Koppens, F. H. L.; Chang, D. E.; Garcia de Abajo, F. J. Nano Lett. 2011, 11, 3370- 
3377. 



15 



(8) Liu, F.; Choi, J. Y.; Seo, T. S. Biosensors and Bioelectronics 2010, 25, 2361-2365. 

(9) Dong, H.; Gao, W.; Yan, F.; Ji, H.; Ju, H. Anal. Chem 2010, 82, 5511-5517. 

(10) Chang, H.; Tang, L.; Wang, Y.; Jiang, J.; Li, J. Anal. Chem 2010, 82, 2341-2346. 

(11) Zhang, M.; Yin, B.-C; Wang, X.-F.; Ye, B.-C. Chemical Communications 2011, 4Z 
2399. 

(12) Chen, Q.; Wei, W.; Lin, J.-M. Biosensors and Bioelectronics 2011, 26, 4497-4502. 

(13) Liu, C; Wang, Z.; Jia, H.; Li, Z. Chemical Communications 2011, 47, 4661. 

(14) Wang, Y.; Li, Z.; Wang, J.; Li, J.; Lin, Y. Trends in Biotechnology 2011, 29, 205-212. 

(15) Nikitin, A. Y.; Guinea, F.; Garcia- Vidal, F.; Martin-Moreno, L. Phys. Rev. B 2011, 

84. 

(16) Swathi, R.; Sebastian, K. J. Chem. Phys. 2008, 129, 054703. 

(17) Gomez-Santos, G.; Stauber, T. Phys. Rev. B 2011, 84, 165438. 

(18) Velizhanin, K. A.; Efimov, A. Phys. Rev. B 2011, 84. 

(19) Stohr, R. J.; Kolesov, R.; Xia, K.; Reuter, R.; Meijer, J.; Logvenov, G.; Wrachtrup, J. 
ACS Nano 2012, 6, 9175-9181. 

(20) Chen, Z.; Berciaud, S.; Nuckolls, C; Heinz, T.; Brus, L. ACS Nano 2010, 4, 2964-2968. 

(21) Treossi, E.; Melucci, M.; Liscio, A.; Gazzano, M.; Samori, P.; Palermo, V. Journal of 
the American Chemical Society 20 09, 131, 15576-15577. 

(22) Kim, J.; Cote, L. J.; Kim, F.; Huang, J. J. Am. Chem. Soc. 2010, 132, 260-267. 

(23) Huang, P.-J. J.; Liu, J. Small 2012, 8, 977-983. 

(24) Heinz, T.; Chen, C; Ricard, D.; Shen, Y. Phys. Rev. Lett. 1982, 48, 478-481. 

16 



(25) Konstantatos, G.; Badioli, M.; Gaudreau, L.; Osmond, J.; Bernechea, M.; de Arquer, F. 
P. G.; Gatti, F.; Koppens, F. H. L. Nat. Nanotech. 2012, 7, 363-368. 

(26) Palik, E. D. Handbook of Optical Constants of Solids; Academic Press: San Diego, CA, 
1998. 

(27) Blanco, L. A.; Garcia de Abajo, F. J. Phys. Rev. B 2004, 69, 205414. 

(28) Barnes, W. L. J. Mod. Opt. 1998, 45, 661-699. 

(29) Lu, C.-H.; Yang, H.-H.; Zhu, C.-L.; Chen, X.; Chen, G.-N. Ang. Chem. 2009, 48, 
4785-4787. 

(30) Glauber, R. J.; Lewenstein, M. Phys. Rev. A 1991, 43, 467-491. 

(31) Wunsch, B.; Stauber, T.; Sols, F.; Guinea, F. New J. Phys. 2006, 8, 318. 

(32) Abramowitz, M.; Stegun, I. A. Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, ninth dover printing, tenth gpo printing ed.; Dover: 
New York, 1965. 



17 



